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ABSTRACT 

Aims. The aim of the project is to characterise both components of the nearest brown dwarf sytem to the Sun, WISE 
J104915.57-531906.1 (also proposed as Luhman 16AB) at optical and near-infrared wavelengths. 

Methods. We obtained high signal-to-noise intermediate-resolution (R ~ 6000-11000) optical (600-1000 nm) and near-infrared 
(1000-2480nm) spectra of each component of Luhman 16AB, the closest brown dwarf binary to the Sun, with the X-Shooter in¬ 
strument on the Very Large Telescope (VLT). 

Results. We classify the primary and secondary of the Luhman 16 system as L6-L7.5 and T0±1, respectively, in agreement with 
previous measurements published in the literature. We present measurements of the lithium pseudo-equivalent widths, which appears 
of similar strength on both components (8.2+1.0A and 8.4+1.5A for the L and T components, respectively). The presence of lithium 
( 7 Li) in both components imply masses below 0.06 M Q while comparison with models suggests lower limits of 0.04 M 0 . The detection 
of lithium in the T component is the first of its kind. Similarly, we assess the strength of other alkali lines (e.g. pseudo-equivalent 
widths of 6-7 A for Rbl and 4-7 A for Csl) present in the optical and near-infrared regions and compare with estimates for L and 
T dwarfs. We also derive effective temperatures and luminosities of each component of the binary: -4.66±0.08 dex and 1305^{|° 
for the L dwarf and -4.68±0.13 dex and 1320^511 for the T dwarf, respectively. Using our radial velocity determinations, the binary 
does not appear to belong to any of the well-known moving group. Our preliminary theoretical analysis of the optical and 7-band 
spectra indicates that the L- and T-type spectra can be reproduced with a single temperature and gravity but different relative chemical 
abundances which impact strongly the spectral energy distribution of L/T transition objects. 
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1. Introduction 


Since the discovery of the first brown dwarfs in 1995 
dNakaiima et alJ 1 1 9951 iRebolo et al.l 1 19951) . the field of sub- 
stellar research has made enormous progress with the dis¬ 
covery of more than 1000 nearby ultracool dwarfs, defined 
as objects with sp ectral types later than M7, which include 
L (~ 1300-2200 K: iMartfnet all 1 1 9991: IKirknatrick efafl 120001: 
iBasri et a ll 120001: iLeggett et al. 12000 ) and T dwarfs (~1300- 
600 K: IBurgasser et alJ 120061) , The coolest b rown dwarfs ever 
found to date, originally nicknamed Y dwarfs (fkirknatric k et alJ 


1 19991) have bee n announced by t he W ide Infrared Survey 
Explorer (WISE; WrighfetaL 2010|) team (ICushing et al.l[201 It 


IKirknatrick et al.l 20121: iTinnev et al1l2012l) . They have temper¬ 

atures estimated to 500-300K and mass es below 0.01 M 0 , ac¬ 
cording to state-of-the-art models (ICushing et al J|20 fTI) . 


During the past decade, e Indi B was the closest brown 
dwar f binary to the Sun , located at 3.626+0 . 009 p c from the 
Sun (IScholz et al.l 120031: iMcCaughrean et aTI [2004). It has a 
mean projected physical separation of 2.65 au (0.75 arcsec) and 


* Based on observations collected at the European Southern 
Observatory, Chile, under DDT programme 290.C-5200(B) (PI Lodieu) 


it is located at ~1 500 au from e Indi A ( Torres et al . 200 6: 
Ivan Leeuw en 2007). It is the best studied p air of brown dw arfs 
with the high est qualit y datas et to date dKing et alJl2010l) . In 
March 2013, iLuhmanl (1201 3l) announced the discovery of a 
nearby brown dwarf with an optical spectral of L8 at a dis¬ 
tance of 2.00+0.15 pc, WISE J104915.57-531906.1A, resolved 
as a close binary (physical separation of 3 au) in the /-band. 
iMamaied 12013l) proposed to call this new nearby brown dwarf 
binary Luhman 16AB due to its proximity to the Sun. The de¬ 
tection of the lithium in absorp tion at 6708A in the optical spec¬ 
trum of the primary ILuhmanl (12013l) and later in both compo¬ 
nents unambiguously places the system in the substellar regime 
(iFahertv et al.ll2014l) . It is brighter than eIndi B by 1.5 mag in J 
and 2 mag in 7. It is the third closest system to the Sun, after the 
Centauri system and Barnard’s star. Thus, it represents the best 
substellar system amenable for detailed characterisation of its 
spectral energy distribution and a unique target to understand the 
chemical processe s at pla y at low temperatures and accross the 
L/T transition (IBurgasser et al.|[2014l) . This system also provides 
a rare opportunity to obtain high-resolution and high signal-to- 
noise spectroscopy at optical and infrared wavelengths to test the 
mass-luminosity-age relation predicted by state-of-the-art mod- 
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Table 1 . Observing logs of the VLT X-shooter observations and information on the wavelength range, resolution, and properties of 
the three arms of the X-shooter spectrograph. 


Arm 

Detector 

Slit 

arcsec 

Resolution 

A range 

jjm 

ExpT 

sec 

Date 

DDMMYY 

Time 

hh:mm 

Airmas 

UVB 

4096x2048 E2V CCD44-82 

O'.'8 

6200 

0.30-0.56 

360 

11062013 

00:31-01:53 

1.29-1.54 

VIS 

4096x2048 MIT/LL CCID20 

O'.'7 

11000 

0.56-1.02 

330 

11062013 

00:31-01:53 

1.29-1.54 

NIR 

2096x2096 Hawaii 2RG 

076 

6200 

1.02-2.48 

6x52 

11062013 

00:31-01:53 

1.29-1.54 


els as well as the chemistry, vertical mixing, abundances of the 
various species present in sub stellar atmo spheres, the role o f 
clouds (ICrossfield et alJ 120141), var iability dGillon et al.l 1201 3l) . 
and search for planets dBoffin et al.ll2014l) . 

In this paper we present medium-resolution optical and 
near-infrared spectra for each component of the nearest brown 
dwarf binary to the Sun, Luh man 16AB (or i ginally known 
as WI SE J104915.57-531906.1: lLuhmanll20 1 3l) . iKniazev et alJ 
(1201 3l) presented resolved optical spectros c opy co vering the 
670-900 nm range, while Bur gasser et alJ (120 1 3h published 
low-resolution (120- 300) near-infrared (70 0-2500 nm) spectra 
for each component. iFahertv et al.l ((20141) presented medium- 
resolution optical (R~4000) and near-infrared (R~8000) of each 
component of the system during the writing of this manuscript. 
In Section [2] we describe the observations conducted in service 
mode with X-shooter on the European Southern Observatory 
(ESO) Very Large Telescope (VLT) and the associated data re¬ 
duction. In Section [3] we assign spectral types to each compo¬ 
nent of the system and discuss the strength of the various al¬ 
kali lines, including lithium. In Section |4] we derive physical 
parameters such as effective and luminosity for each component 
of the binary. In Section [3 we derive the radial velocity of the 
system and discuss its membership to nearby moving groups. 
Synthetical spectra of the alkalis are computed and compared to 
the Xshooter data in Section[6] 


2. Optical to infrared spectroscopy 

We carried out spectroscopy from the UV- to the K- 
band with the X-shoote r cross-dispersed ec helle spectrograph 
(iD’Odorico et al.l 120061: IVernet et al.l 1201 ll) mounted on the 
Cassegrain focus of the Very Large Telescope (VLT) Unit 2. 
Observations were carried out in service mode on 11 June 2013 
under clear skies and a seeing between 0.6 and 0.85 arcsec to 
ensure that both components of the binary system were resolved 
over the full wavelength range covered by X-Shooter. More de¬ 
tails on the properties of the X-shooter spectrograph and the logs 
of the observations are provided in Table Q] 

The acquisition of the target was made through an optical 
CCD with the Sloan z-band filter, which is not problematic con¬ 
sidering the brightness of the source. We set the individual on- 
source integration times to one expossure of 360 sec, one of 330 
sec, and six of 52 sec in the UVB, VIS, and NIR arms, respec¬ 
tively, repeated five times with a maximum nod-throw of 5" to 
correct for the sky contribution (mainly) in the near-infrared. The 
field-of-view and the slit were both rotated by a position angle 
of -47 degrees to place both components in the slit. Hence, the 
observations were not conducted with the slit orientated along 
the parallactic angle. We caution the use the optical spectrum of 
each component of the binary for accurate spectral classification 
because of the effect of differential refraction particularly at high 
air masses. However, the relative fluxes between components re¬ 
main unaffected because both sources are separated by about 1'.'5 


and have similar spectral types, hence, they are supposedly af¬ 
fected the same way. 

Data reduction was performed for the VIS and NIR arms 
separately. For the present paper, we discarded the UVB arm 
because no obvious signal is registered from any of the com¬ 
ponents. We applied standard procedures for optical and vis¬ 
ible wavelengths, which include sky subtraction and flat-field 
correction. All steps were carried out with routines within the 
IRAFL 1 ] (Image Reduction and Analysis Facility) environment 


_ ma gi 

Todvll 1986111993b . Spectra were optimally extracted with the 
APALl command of the ECHELLE package in IRAF. Each 
pair member was extracted separately using the same aperture 
size to define the extraction and background regions on the de¬ 
tectors. A third-order polynomial was employed per echelle or¬ 
der to trace the spectra across the detectors. Wavelength cali¬ 


bration was performed to a typical accuracy of +0.05 A using 
ThAr (VIS) and ArKrXeNe (NIR) arcs taken with the same in¬ 
strumental configuration as our targets. The spectra of the B9V 
star HIP 56470 observed immediately after Luhman 16AB with 
an air mass within +0.1 were used for division into the corre¬ 
sponding VIS and NIR science spectra. We previously removed 
the hydrogen lines intrinsic to this hot star. Finally, we multi¬ 
plied the science spectra by a blackbody of temperature 10700 
K corresponding to B9V. Individual echelle orders were merged 
by matching the counts in common wavelength ranges from one 
order to the next. We observe a peak and a dip around 2.08- 
2.09 microns in the spectra of the components of the system not 
seen in other T dwarfs, as a potential artefact in our X-shooter 
data. Though we could not remove it, it does not affect our sci¬ 
entific interpretation, as a potential artefacts. However, it does 
not affect our interpretation. The other regions of the spectra is 
not affected by any other strong artefact(s). The final optical and 
near-infrared one-dimensional spectra are displayed in Figures [I] 
and [2] The latter al so includes known L and T dwarf templates 


from the literature (Knaou et al.l 2004: 

Golimowski et al. 2004a: 

Burgasser et al. 2004; Chiu et al. 2006 

Looper et al][2007]). 


3. Spectroscopic properties of each component 

3. 1. Spectral types 

In Fig. Q] we compare the full spectra (660-2500 nm), nor¬ 
malised at ~1300nm, of the L and T components (Luhman 16A 
and Luhman 16B, respectively). We see that the L component 
is relatively brighter than the T component bluewards of 1000 
nm. The overall fluxes are very similar in the /-band while the 
L component dominates in the H and K bands. We confirm that 
Luhman 16A is of earlier type than Luhman 16B. 

We assign spectral types to each component of the 
binary, using direct comparison with known L and T 

1 IRAF is distributed by National Optical Astronomy Observatories, 
whcih is operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foundation. 
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Fig. 1. Comparison of the X-shooter optical and near-infrared spectra of the L (black line) and T (red line) components of the 
Luhman 16AB brown dwarf binary. Both spectra are normalised to unity at 1300 nm to highlight the differences in spectroscopic 
features between both components. Both spectra are also corrected for telluric contribution. The main spectral features present in L 
and T dwarfs are highlighted. 



spectral templates both at optical and near-infrared wave¬ 
lengths availa ble at the S peX a rc hival and Sandy L eggett’s 
homepag jl d Chiu et alJ 120061: iGolimowski et al.1 l2004at 
Knan n et al.l l2004ll . We classified each component sep¬ 
arately in the optical and the near-infrared. The optical 
spectrum of the T component is b est fit by the L 9 tem¬ 


plate 2MASSJ03105 9 86+1648155 dkirknatrick et a! 


i Bureasser et al.l 20061 iReid et all 1200 It iLeggettetaT 


2000 ; 


2002 ; 


Geballe et al. 2002) whereas its near-infrared is best repro- 


duced by th e TO te mplate SPSS J 04234 8 57-0414035 (TO; 


Leggett et al.l l2000t ICieballe etal.1 l200l ICruz et al.l [2003; 


Burg asser et al.l120*061) with an uncertainty of one subtype (Fig. 


3. The L component is not well reproduced by any of the L 
dwarf spectral templates over the full spectral energy range. 

Our final spectral classifications, T0+1 for Luhman 16B and 
L6-L7.5 to Luhman 16A (Fig. [2}. agree with previous stud¬ 
ies, such aslBurg asser et al] (1201 3T L7.5+0.5 and TO.5+0.5) and 
I Kniazev et al.l (I2013L L8+1 and T1 ±2). 


3.2. Molecules 

In the optical spectra of both components, we observe a wide 
pressure-broadened potassium doublet around 760 nm. We can 
also see the FeH, CrH, and FLO absorption bands at ~860nm, 
~870 nm, and ~930 nm, respectively, pointing towards cool tem- 


2 http://pono.ucsd.edu/~adam/browndwarfs/spexprism/ 

3 http://staff.gemini.edu/ sleggett/LTdata.html 


peratures and la te spectral types (Kirkpatri ck et alJ|l999L [2000: 
IReid et alJl2000h . 

In the near-infrared, we observe a stronger methane absorp¬ 
tion band at around 1.1-1.2 microns in the T component, in¬ 
dicating a potentially cooler temperature. By comparison, the 
shape of the spectra in the H -band appear very similar, point¬ 
ing towards a small difference in spectral type. We detect the 
onset of methane (CH 4 ) absorption at 2.12 /rm in Luhman 16B, 
which clearly indicates a T spectral type. Additionally, we find 
CO in absorption at 2.3 fim in both L and T components. The co¬ 
existence of CH 4 and CO hints at a peculiar carbon chemistry. 

3.3. Lithium 7 C Li) 

The presence of lithium in the atmospheres of fully convec¬ 
tive dwarfs cooler than M 6 has been widely used as a solid 
criterion for substellarity dRebolo et al .11 19921 : iBasri et afll2000b 
Kirkpatri ck et alJ[20081 and references therein). We detect the 
lithium resonance doublet at 670.782 nm in both components, 
thus confirming their brown dwarf nature and masses <0.060 M e 
independently of the age of the binary. This is the first detection 
of Li 1 absorption in a T dwarf (also pointed out bv lFahertv et al.l 
120141 ')'). We measured the equivalent width of the spectral fea¬ 
ture (note that the doublet is not resolved) with respect to the 
objects’ relative continuum or pseudo-continuum modulated by 
the strong absorptions due to Nai to the blue and Ki to the red. 
We therefore refer to pseudo-equivalent widths (pEWs). The in¬ 
tegration of the line profile is always performed over the range 
669.45-672.25 nm. Our measurements and their associated un- 
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Fig. 2. Top panels: VLT X-shooter optical (left) and near-infrared (right) spectra corrected for telluric lines for the L dwarf com¬ 
ponent of Luhman 16AB. Bottom panels: Same but for the T component of Luhman 16AB. Overplotted as red lines are spectral 
standards downloaded fr om the SpeX arch ive an d Sandy Leggett’s w ebpage: 2MASS J15074769-1627 38.6(L5-L5.5; RekLetaLl 
2000: iDahn elaD 12002: Kirknatr icketaP l2000t iKnapp et allhool . 2MASSsJ0850359+105 716 (L6: lKirkpatricketal.il 19m 


2MA SSI JO 103320+193536 (L7: iKirknatric kTt alJl2000h, 2MASS J 16322911 + 1904407 (L7.5: Eeggett et al 


2002) . 2M ASS .103 105986+1648155 (L9: [Kirknatri dTet al.ll2000t iBurgasser et all l2006t iReid et al.l 1200 


Geba lle et al.1 12002). SPSS J04234857- 0414035 (TO: iLeggettetal 


2006), SPSS.T08371721- 0000180 fTl: iLeggett et alJl2000l 1 20021 


J125453.90-012247. 4(T2: iLeggett et all20Q 
860 nm and 1300 nm to ease the comparison. 


200(J Geballe et al. 1120021: ICruz et al. 20031: Burgasser et al.1 


2002 : 

Legge 


Burgasser et al.1 2006T Kirkpatrick et alJ 2' 


Geballe et al. 

gett et al. 2002: 


and SPSSp 

Burgasser et all2004t) . All optical and near-infrared spectra are normalised at ~840- 


certainties (determined as the scatter of the measurements ob¬ 
tained by placing the pseudo-continuum up and down within the 
photon noise of the observed spectra) are provided in Table [2] 
Interestingly, we find that both the primary and the sec¬ 
ondary have similar lithium absorption intensities of pEW ~ 8- 
9A within the error bars despite their differing spectral types. 
Figure [3] illustrates the Li i profiles of Luhman 16A and Luhman 
16B. This contrasts with the Li i pEWs reported bv lFahertv et al.l 
(2014). These authors used spectra of resolution R ~ 4000 at I- 
band taken 1.5 months earlier than our X-shooter data and de¬ 
rived that the T-type component showed a lithium spectral sig¬ 
nature smaller (pEW = 3.8 + 0.4A) than that of the L-type com¬ 
ponent (pEW = 8.0 ± 0.4A). Whereas our pEW measurement 
agrees with iFahertv et al.l (2 0141 ) result at the 1-cr level for the 
L dwarf, there is a discrepancy by a factor of ~2 for the T 
component. We repeated our procedure on the optical spectra 
kindly provided by Jackie Faherty and measured similar pEWs 


in Luhman 16A and Luhman 16B of the order of 7-8 A, in ex¬ 
cellent agreement with our X-shooter pEWs. Thus, the differ¬ 
ent pEWs determinations are likely due to distinct measurement 
techniques rather than to variability. 


Our Li i pEWs and those of IFahertv et alJ ( 2014 ) are de¬ 
picted as a function of spectral type in Fig. [4] To put our 
values in context, we also include pEWs of field late - M, L, 
and T dwarf s given bvlTinnevI (11998 ), IKirknatrick et alJ (1 19991 
2000|^ 2008 k_ Pavlenko igtakl (l2007£ ), ICruz et al.l (l2()09l) . and 


Zapatero Osorio et al 


Imi . As shown in Fig. |4j the mea¬ 
surement of Luhman 16A is consistent with the pEWs of L7- 
L8 dwarfs with lithium detections. There is a large scatter in 
Li i pEWs of L dwarfs possibly due to the uncertainties asso¬ 
ciated with the measurements and spectral types, the unequal 
spectral resolutions and ways of deriving pEWs by the various 
groups, distinct states of lithi um depletion, variability, and/o r di¬ 
verse metallicity. However. IKirknatrick et al.l (1200(1 [2008) dis- 
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Table 2. Pseudo-equivalent widths (in A) of atomic features. 


Line" 

(nm) 

Luhman 16A 

Luhman 16B 

A4* 

(nm) 

Li 14670.782 

8.2+1.0 

8.4+1.5 

669.45-673.50 

Rb i 4780.023 

6.2+0.7 

6.7+0.8 

779.04-781.44 

Rb i 4794.760 

6.7+0.7 

7.1+0.7 

793.12-796.71 

Na 1 4818.326 

1.0+0.2 

0.5+0.1 

818.00-819.00 

Nai 4819.482 

1.3+0.2 

0.8+0.1 

819.00-820.15 

Cs i 4852.115 

6.5+0.7 

7.4+0.7 

849.85-854.44 

Csi 4894.348 c 

4.2+0.5 

5.0+0.5 

893.65-895.07 

Ki 41169.022^ 

5.3+0.5 

7.5+0.5 

1167.93-1170.59 

Ki 41173.284'' 

5.4+0.5 

7.1+0.5 

1176.34-1178.91 

Ki 41243.227 

3.0+0.5 

4.6+0.5 

1241.15-1245.00 

Ki 41252.213 

4.0+0.5 

6.1+0.5 

1250.25-1254.25 


Notes. (a) Catalogue wavelengths are given in the air system. 
(i) Wavelength interval over which the line profile is integrated. 
(c) Blended with FeH absorption . {d) Blended with H 2 0 absorption. 



Fig. 3. Spectral profiles of the Lii resonance doublet at 670.782 
nm seen in absorption. The lines of the two objects are similar 
in strength. The relative fluxes of both pair members are normal¬ 
ized to the mean pseudo-continuum around the Li i signature for 
a proper comparison. The data are free of telluric contribution. 
No smoothing has been applied. 


cussed that the strength of the Lii 670.782 nm line increases 
with later spectral type, peaks at L5-L6, and declines toward 
cooler types. As pr edicted by the theory of stellar and substel- 
lar chemistry (e.g., iBurrows & Shard 1199 9: [Lodders & Feglevl 
20061 ) . lithium is in atomic form above ~ 1500-1700 K at 1 
bar pressure but is converted into molecules (LiOH, LiCl, LiH, 
...) at lower temperatures. The resulting effect on the ob¬ 
served spectral energy distribution is the depletion of the atomic 
lithium abundance and the decreasing intensity of the Lii ab¬ 
sorption at 670.782 nm toward very low temperatures and late 
spectral types. Nevertheless, this feature is still present and 
remains detectable with relatively high strengths in early-Ts, 
as demonstrated by our observations of Luhman 16B and in 



Fig. 4. Li 

T dwarfs. Data 


pseudo-equivalent widths (pEWs) of late-M, L, and 
tak en from iTinnevI (1 19981) Kirkpatrick et alJ 

. IPavlenko et alJ (l2007a). Cr uz et alT (2009). 


and lZanatero Osorio et al.l(l2014l) are shown with black dots (de 


tections) and arrows (upper limits). Our pEWs of Luhman 16A 
(L6.75) and Luhman 16B (TO) are indicated with the red dots, 
and lFahertv et alJ (1201 4lf measurements are plotted as blue solid 
triangles. Values from lFahertv et alJ (1201 4l) for Luhman 16A are 
slightly shifted in spectral type and pEW by +0.1 subtype and 
+0.1 A, respectively, for clarity. 


agreement with the predi ctions made by IPavle nko et al.l (|2000), 
lAllard et~afl d200ll) . and IBurrows et al.l d2002ah . Furthermore, 
the pEW = 8.4+ 1.5A of Luhman 16B is compatible with the 
computed intensity of the Li i resonance doublet for 77 n = 1000- 
1200 K, full lithium preservation, and moderate dust opacity 
given in Table 5 of lPavlenko et al.l d2000l) . 

Among the Lii subordinate lines, the one at 812.64 nm 
would be more easily detectable because there is more flux (bet¬ 
ter signal-to-noise data) in our targets than at 610.36 nm. The 
detection of one of these lines would provide a direct evidence 
of the full preservation of lithium in the dwarfs’ atmospheres, 
since the formation of the subordinate lines is extremely sen¬ 
sitive to the atomic abundance. Unfortunately, the detection of 
these lines in absorption also depends on the presence of dust 
opacity, which considerably reduces the line strength. We can 
impose an upper limit of 30 mA on the Li 1 4812.64 nm pEW for 
both components. A detailed spectral analysis of the binary data 
would be required to confirm whether lithium is fully preserved 
or depleted in both components. 

3.4. Lithium 6 ('Li) 

Altough the high superficial gravity in brown dwarfs produces a 
large broadening of the alkali lines, in particular of the 7 Li reso¬ 
nance doublet, it is worth to study the potential presence of the 
more fragile 6 Li isotope in low- mass brown dwarfs. According 
to theoretical computations (e.g. lNelson et al.’f 19931) . 6 Li, which 
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Fig. 5. Spectral profiles of Nai /1818.326 and /1819.482 nm 
subordinate lines. The relative fluxes of the two objects have 
been normalized to the mean pseudo-continuum redwards of 
T819.482 nm. Spectra are telluric corrected and no smoothing 
is applied. 

is burnt at lower temperatures than 7 Li, would be fully preserved 
in brown dwarfs of mass lower than 40 Mj up . Given the poten¬ 
tial masses and the brightness of the Luhman 16AB binary, these <u 
brown dwarfs offer the interesting opportunity to investigate the 3 
presence of 6 Li. The isotopic shift between the doublets of 7 Li % 
and 6 Li is only of 0.16A and therefore requires a fine spectral 
resolution and high signal-to-noise. 

A detection of 6 Li would clearly hint to a mass below 40 
Mj up . The lithium isotopic ratio in meteorites is 7 Li/ 6 Li ~ 12, 
and es tima tes in the interstellar medium give similar values 
iLemoine et all 4993 ), so the presence of the 6 Li isotope is ex¬ 
pected to introduce a weak asymmetry in the red wing of the 
lithium doublet. We investigated if there is any evidence for 
such asymmetry in the spectrum of Luhman 16A (with a bet¬ 
ter signal-to-noise than Luhman 16B) by subtracting a Voigt and 
a Gaussian function from the observed profile. We found flat 
subtractions with no significant deviation (particularly in the red 
wing of the absorption) at the level of 1/8 of the Li relative flux 
(1 -cr) that could be attributed to the presence of 6 Li. Higher spec¬ 
tral resolution and improved signal-to-noise spectroscopy of the 
lithium feature is required to impose significant constraints on 
the presence of the lightest lithium isotope and in turn on the 
masses of these brown dwarfs. 

3.5. Other alkali lines 

Rb i at 780.023 and 794.760 nm, and Cs i at 852.115 and 894.348 
nm are clearly seen in the X-shooter optical spectra with pEWs 
given in Table[2] We measured these pEWs and their error bars in 
the same manner as for lithium. Fig. [6] displays an enlargement 
of the X-shooter spectra around the Rb i and Cs i lines; we do 
not smooth the data for this comparison. The intensity of the 
Rb i atomic signatures is quite similar within the error bars for 
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Fig. 6. Spectral profiles of Rbi (top) and Csi (bottom) atomic 
lines of Luhman 16AB. The relative fluxes of the two objects 
have been normalized to the mean pseudo-continuum around 
Rbi and the Csi 4852.11 nm lines for a proper comparison. 
These data are corrected for telluric contribution. The effect of 
the different T e g between the components becomes apparent in 
the slope of the pseudo-continua in the bottom panel. The Cs i 
4894.35 nm line is blended with the absorption due to CrH and 
FeH. 


both the L and T components (see top panel of Fig.[6]»; however, 
Cs i appears slightly stronger in the T dwarf (bottom panel of 
Fig. [6]), which may be a signpost of its cooler temperature. 

Contrary to the Lii resonance line (Section [33]). our pEWs 
of Rb i and Cs i agree with the equivalent widths published 
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Fig. 7. Csi /i852.115 nm pseudo-equivalent widths (pEWs) 
of L and T dwarfs. Our pEWs of Luhmanl6A (L7.5) 
and Luhman 16B (TO) are indicated with the red dots, 
and iFahertv et alJ (l2014l> measu rements are plot t ed as blue 
solid triangles. Data taken from Kirkpatri ck et alJ (l2000h and 
Burgasse r et alJ (120031) are depicted with black solid circles, and 
our measurements over optical spectra available in the literature 
are plotted as open circles. 


bv IFahertv et al.1 J2014i) within l-cr the uncertainty. The mea¬ 
sured pEWs are also consistent wit h the values found for 
late-Ls and early-Ts in the literature (Kirkpatrick et al. 2000; 
(Burgasser et all 120031) . as it is illustrated in Fig. [7] where the 
strength of Cs i /1852.115 nm is shown as a function of spectral 
type. We choose this Cs i feature because it is located far red 
from the strong K i absorption in the optical and it is not blended 
with other atomic or molecular absorptions. To complete the dis¬ 
played data, we measured pEWs over optical spectra of field 
L6-T6 dwarf s available to us and publ is hed bvlKirk patrick e t alJ 
(19991 20001), IBurgasser et all (120001) . iReidetal* (12001 ). and 
ICruz et alJ (120031) . These measurements are plotted as open cir¬ 
cles in Fig. |7| Similarly to what is observed for lithium, the in¬ 
tensity of the cesium line rapidly increases toward later spec¬ 
tral type and reaches a maximum between TO and T5, this is 
at temperatures cooler than the lithium case. The marked inten¬ 
sity decrease at >T5 is likely due to the conversion of atomic 
Cs i into molecules as a natural result of the cool atmospheres’ 
chemistry. This property is in good agreement with the the¬ 
ory, which predicts that rubidium and cesium disappears in 
monatomic form at sufficiently low temperatures that are 200- 
300 K smaller than the one f or lithium (e.g.. lBurrows & Sharp 
fl999t iLodderiT&F eglev 20061). 


We detect Nai T818.326 and 4819.482 nm subordinate lines 
in the spectra of Luhman 16A and Luhman 16B. The detection 
is highlighted in Fig.0 We report their pEWs in Table [2] Both 
components have Nai p EWs greater than 0.5A, which contrasts 
with the claim bv iFahertv et al.l (i2014 ). These authors imposed 
an upper limit on the strength of Nai absorption of <0.5A. This 
discrepancy can be explained by the different spectral resolu¬ 


tion of the data and the fact that iFahertv et al. ( 2014 ) spec¬ 
tra are not corrected for fringing and telluric contribution (im¬ 
portant at these wavelengths), as acknowledged by the authors. 
Additionally, the optical spectra of lKniazev et al.l ( 2013 ) are not 
free of telluric absorption; therefore, these authors were not able 
to detect this feature despite the good quality of their data. The L 
component shows stronger Nai lines than the T member of the 
pair, probably because the pseudo-continuum of the latter ob¬ 
ject is more depressed due to a lower temperature and a more 
vigorous red wing of the Ki doublet centered on 768.2 nm. 
Interestingly, the presence of Nai at 1138.145 and 1140.378 
nm in the near-infrared is not obvious in either pair member, 
in agreement with iFahertv et al. ( 201 4). Other early T dwarfs, 
like e Indi Ba, do not show these lines in their high quality 
spectra (King et al.ll2010h . Thus, we do not confirm the ne ar- 
infrared detection of Nai reported by Burgasser et al.l (1201 3l) . If 
any, only the reddest line of the two (1140.378 nm) may display 
a marginally resolved signature with pEWs = 1.3A (L) and 1.5A 
(T), and no trace of the blue line. Note that this region is severely 
affected by HiO absorption. 

Potassium is detected at optical and near-infrared wave¬ 
lengths. The resonance doublet at 766.4911 and 769.8974 nm 
extends several 1000A (typical of late-L and T dwarfs) and im¬ 
prints the overall shape of the visible spectra. We report the 
pEWs of the near-infrared Ki lines in Table [2] The strength of 
the optical and near-infrared Ki lines (766.4911 and 769.8974 
nm, 1.169, 1.173, 1.243, and .1252 yum) is systematically larger 
in Luhman 16B than in the L dwarf, which may suggest a higher 
potassium abundance, a lower temperature, depletion of potas¬ 
sium in condensates, and/or a pseudo-continuum affected by less 
opacity at near-infrared wavelengths. The properties of the near- 
infrared K i lines are widely discussed in comparison with the 
spectra of other field L and T dwarfs bv lFahertv et al.l (|2014|) . 

It is known that the alkali lines have a profound correlation 
with gravity (atmospheric pressure) and metallicity at the low 
temperatures of the binary: the higher the gravity and the atomic 
abundance, the stronger the spectral features. With the only ex¬ 
ception of lithium, the similarity of the strength of the alkali lines 
between Luhman 16A and B and field ultra-cool dwarfs of re¬ 
lated spectral types indicates that the components do not deviate 
significantly in age and metallicity with re spec t to the field, in 
agre ement with the age derivations by Burgas ser et al.1 (1201 3l) 
and lFahertv et alJ (120141) . 


4. Effective temperatures and luminosities 


We determined the luminosity of the L and T c omp onents using 
,/-, H - , and K-band bolomet r ic corrections (Golimowski et al.l 
l2004bt iNakaiima et al.l l2004t IVrba et aTI 120041) valid for field 
dwarfs and the spectral intervals L6-L8 (Luhman 16A) and L9- 
T1 (Luhman 16B), which accounts for the spectral types of the 
pair members and their associated uncertainties. Note that the 
H -band bolometric correction (BC) shows the smallest disper¬ 
sion among the three bands. We also used a solar bolomet¬ 
ric luminosity of 4.73 mag and the trigonometric distance of 
2.02 ± 0.02 pc determined for the system bv lBoffin et al.l ( |2014|) . 
The near-infrared m agnitude s of each pair mem ber provided by 
Burg asser et al.l (120131) and lKniazev et al.l (l2013h were averaged. 
We derived the following values: log L/L e = -4.68 ± 0.08 (L 
component) and -4.66 ± 0.08 dex (T component). The luminos¬ 
ity error bars consider the uncertainties in the photometry (typ¬ 
ically +0.03 mag), the distance modulus (±0.02 mag), and the 
corresponding BCs (typically ±0.15 mag for J and K , and ±0.08 
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Table 3. Masses, radii, and effective temperatures (T e ff) for each component of Luhman 16 for different ages along with the error 
bars from the interpolation of the models. 

Note: a this age yields a mass that is not compatible with lithium preservation 


Age 

Gyr 

Mass 

M 0 

Luhman 16A 
Radius 

Rg 

T e ff 

K 

Mass 

M 0 

Luhman 16B 
Radius 

Rg 

T e ff 

K 

0.5 

0.028+0.003 

0.102+0.005 

1223+60 

0.029+0.003 

0.102+0.005 

1238+60 

1 

0.041+0.003 

0.093+0.005 

1280+60 

0.041+0.003 

0.093+0.005 

1296+60 

2 

0.055+0.003 

0.086+0.005 

1328+60 

0.055+0.003 

0.086+0.005 

1344+60 

5“ 

0.067+0.003 

0.082+0.005 

1368+60 

0.067+0.003 

0.082+0.005 

1383+60 


for H). Our determination for the L component broadly agrees 
with Fahertv et al.l (120141) . who found log L/L 0 = -4.67 + 0.04 
dex after the integration over their optical (0.6-0.9 /urn) and near- 
infrared (0.8-2.5 fm t) spectra supplemented with synthetic spec¬ 
tra for wavelengths longer than 2.5 /mi. This supports the relia¬ 
bility of the BCs for these spectral types. Regarding the T com¬ 
ponent, our luminosity determination is consistent with that of 
IFahertv et al.l (120141 -4.71 +0.10 dex) at the 1 -cr level, although 
our measurement suggests a slightly brighter luminosity. In a 
similar manner as described in IFahertv et al.l d2014l) . we also in¬ 
tegrated our X-shooter spectra and found that the L and T dwarfs 
have nearly identical luminosity, with the T object being 0.01 
dex more luminous than the L dwarf, in agreement with the lu¬ 
minosities inferred from the BCs. The T component has been 
reported to show larger variability than the L source at optical 
and near-infrared wavelengths with peak-to-peak amplitudes of 
11-13.5 % in continuous observations and s trong night- t o-nigh t 
evol ution according to iGillon et al.l (1201 3l) . iBiller et aD (1201 3l) . 
and iBurgasser et al.l (120141) . This might impact the luminosity 
determination for this particular dwarf. To account for the ob¬ 
served variability amplitude, an uncertainty of ±0.05 dex should 
be added to the error bar in the luminosity derivation of the T pair 
member, thus yielding log L/L 0 = -4.66+0.13 dex. One striking 
result is that, despite the differing spectral classifications, both 
objects, Luhman 16A and B, have consistent L/L 0 values at the 
1-cr level. 

There are several methods to determine the effective tem¬ 
peratures (r eff ) of Luhman 16A and B. One is fitting observed 
spectra using theoretical model atmospheres, as it was done by 
IFahertv et aid d2014t> . 

These authors found a wide range of possible temperatures 
(900-1800 K) for both components, and that diverse models 
yield temperatures that can differ by 500 K for the same source 
and observed data. Another way is to use T e g calibrations as a 
function of spectral type or colour ava ilable in the literature (e.g. 
IVrba et akll2004t fStenhens et al.ll2009l) . Using these two calibra¬ 
tions valid for field dwarfs, we inferred T e « = 1460 ± 100 K for 
the L component, and 1335 ± 100 K for the T component. The 
uncertainties stand for the dispersion of the calibrations as in¬ 
dicated by the authors. A third method consists of applying the 
astrophysical equation 


log L/L 0 = 2 log (R/R 0 ) + 4 log (Teff/Tefr,©), d) 

where r e ff >0 = 5777 K and R/R 0 indicates the object’s radius in 
solar units. This way has the advantage that substellar objects of 
ages greater than ~0.5 Gyr change little in size acco rding to the 
theory of substellar evolution (see review bv IChabrier & Baraffe 
120001) . Particularly, for ages >1 Gyr and brown dwarfs with 
masses >0.02 M 0 , models predict R/R 0 to be in the range 0.076- 
0.102. Considering a mean size of 0.09 R 0 and the bolometric 



Fig. 8. Luminosity versus effective temperature diagram for the 
two components of the Luhman 16AB pair (red open circles). 
Isochrones are plotted as solid lines, and the tracks of 0.020 and 
0.060 M 0 are displayed with dot-dashed and dashed lines. The 
models by Bar affe et al. l d2003l ) are used. Blue triangles stand for 
IFahertv et al.l < 20 14l) determinations. 


luminosities given above, we determined 7 0 r = 1305 + _|Ifi for the 
L component and 1320*^5 for the T component. The generous 
error bars include the uncertainty in luminosity and all possible 
radii between 0.076 and 0.102 R 0 . These error bars would de¬ 
crease by about 100K if we would only consider the uncertain¬ 
ties in size. The determination of the relative luminosity from 
the integration of the spectra has a relative error estimated at 
+/-0.04 dex, which translates into a relative error on the 7 0 r 
of +2.3% or ±30 K at 1300K. Hence, both objects have simi¬ 
lar luminosity within +/-0.04 dex and similar T e g within 30 K 
(lcr). We do not expect different masses or sizes from the the¬ 
ory if the two objects are coeval, which we assumed here. We 
list in Table [3] the masses, radii, and 7 0 r for each component 
of the Luhman 16 system, assuming four different ages (0.5, 
1, 2, and 5 Gyr) along with the error bars from the interpola¬ 
tions of models for each age. We note that the oldest age yields 
a mass that is not compatible with lithium preservation as al¬ 
ready discussed by IBurgasser et all (120141) . Equation Q] and the 
literature T e g —spectral type calibrations provide similar tem- 
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peratures for Luhman 16B, while Eq. Qjyields a cooler temper¬ 
ature for Luhman 16A, which may imply a later spec tral classi- 
fication (e.g. L8) as the one quoted in by (lLuhmanll2013l) and 
iKniazev et al.1 (1201 3l) . Independently of the method used, we 
found that the two pair members have consistent temperatures 
and luminosities, which indicates similar masses and radii. 

We show the bolometric luminosity versus T c n diagram for 
the two components in Figure [8] The temperatures derived using 
Eq.Q]are employed. We overplotted several isochrones (lines in 
colours) from the COND models of lBaraffe et alJ (120031) . rang¬ 
ing from 10 Myr up to 10 Gyr. We also plotted isomasses for 
0.020 and 0.060 M 0 . As it is seen from the Figure, the binary 
appears to be older than 120 Myr with a likely age close to 1 Gyr 
and a mass of ~ 0.040 M e for each pair member. These determi¬ 
nations fully agree with the presence of lithium in the objects’ 
atmospheres ( see Sectionl3.3l). and wi th th e age and mas s range s 
determined by Burgasser et al] d2013h and lFahertv et al.l ([2014). 


5. System velocities and relative radial velocity 


We measured the central wavelengths of the various atomic lines 
indicated in Table[2]to estimate the radial velocity of each com¬ 
ponent. The values of the observed central wavelengths of Rb i, 
Na i, Cs i, and near-infrared K i lines are compared to their nom¬ 
inal air wavelength]]]. The average wavelength shift lies in the 
0.10-0.15 nm range, implying a mean observed velocity of about 
46 kms -1 with a most probable range of 36-56 kms . The 
measurements for both components agree within 10%. We also 
cross-correlated the X-shooter spectra of the L and T members 
against the 1400-K BT-Settl model spectrum computed for so¬ 
lar metallicity and surface gravity of log g = 4.5 (cm s 2 ) by 
1 Allard et~akl ( 2012 ). We used the wavelength interval 775-1350 
nm, which excludes regions of strong telluric absorption and 
poor signal-to-noise ratio and includes atomic and molecular ab¬ 
sorption features due to Rb i, Na i, Cs i, Ki, IGO, CrH, and FeH. 
The resulting observed velocities are 40.9 + 7.6 and 35.1 ±7.2 
kms -1 for the L and T components, respectively. We applied 
a correction of -15.8 km/s corresponding to the diurnal, lunar, 
and annual velocities during the observations on 11 June 2013 at 
UT = 00 h 38 m to obtain the following heliocentric radial veloci¬ 
ties: Vh = 25.1 ±7.6 and 19.3 ± 7.2 kms -1 for Luhman 16A and 
Luhmanl6B. Our measurements and the heliocentric velocities 
(23.1 and 19.5 kms -1 ) reported by IKniazev et alJ d2013l) agree 
within 1 -cr the uncertainty. 

Using our heliocentric radial velocities (mean value of 
22.2 + 7.5 kms -1 for the system), the recent determinations of 
parallax and proper motion b vlBoffin et al.l d2014i ), and the equa¬ 
tions of Johnson & SoderblomF i 19871) . we derive the following 
Galactic space velocities: U = -18.3+2.0, V- -28.7+7.2, and 
W = -6.8+0.7 kms -1 . The uncertainties associated to all three 
Galactic velocities come from the proper motion, parallax, and 
radial velocity errors. The space velocities of Luhman 16AB 
are statistically consistent with the Galaxy young disk kinemat- 
i cs acc ording to the classification made bv lEggei] (1 990 ) and 
iLeggettl ( 1 9921), a result that is in agr eement with the age estimate 
made by Burgasser et alJ (1201 3l) and lFahertv et alJ (1201 4l) and the 
value reported in the previous Section. The UVW velocities do 
not overlap with any close (<100 pc), young ^600 Myr) stel¬ 
lar mo ving group listed i n the catalogues of Zuc kerman & Song 
(120041) andjTorres^talJd2008|), which supports the conclusions 
made bv lKniazev et al.l (l2013l) . 


http://www.nist.gov/pml/data/asd.cfm 


We also measured the relative velocity of the T component 
versus the L component by means of the cross-correlation tech¬ 
nique. We cross-correlated the spectra with and without the tel¬ 
luric contribution. The regions dominated by telluric absorption 
were used to establish the velocity shift corrections to be applied 
for bringing the data of the two components to the same system 
of reference and obtaining an accurate relative radial velocity. 
We derived relative velocities of -1.06 ± 0.50 and -1.42 ± 0.70 
kms -1 from the optical and near-infrared parts of the X-shooter 
spectra, respectively. The associated error bars are smaller than 
those of the heliocentric radial velocities as expected for rela¬ 
tive velocities. The measured mean relative velocity is consis¬ 
tent with the difference between the heliocentric velocities given 
above at the 1-cr level. Our value is also compatible with that 
(-2.5 ± 1.9 kms -1 ) of Kniazev et al.l (1201 3l) obtained from data 
that were acquired 90 days prior to the X-shooter spectra. The 
non-zero relative velocity on two different occasions indicates 
that the binary orbit is likely tilted with respect to the line of 
sight. Precise radial velocity monitoring of the system is required 
to search for unresolved companions (including planets) around 
any component (iBoffin et al,l2014l) . to characterize the binary or¬ 
bit (with an estimated orbital period greater than 15 years), and 
to determine dynamical masses of the components, which will 
in turn become useful for comparison with substellar models. 


6. Comparison with theoretical models 

We compared the observed spectrum of each component of 
the Luhman 16AB system with the BT-Settl synthetic models 
(iAllard et al.l l2012h . We were unable to simultaneously repro¬ 
duce the X-shooter data of the two objects with th e sam e tem¬ 
perature, gravity and metallicity (see also iFahertv et al.l 12014 1). 
In Fig. |9] the BT-Settl spectrum computed for solar abun¬ 
dance, log g- 5.0 and T e s = 1400 K (close to the values de¬ 
rived/discussed in previous sections) is shown in comparison 
with the X-shooter of Luhman 16A and B. However, accord¬ 
ing to IFahertv et alJ ( 120141) , the BT-Settl models can repro¬ 
duce each component of Luhman 16 using different tempera¬ 
tures and surface gravities (i.e., ages), hypothesis unlikely for 
a nearby co-moving, lithium bearing, low-mass binary system. 
Using our code (see next), the fitting of the X-shooter spectra 
yields T e g = 1500 K and log g - 5.0 dex for the L component, 
and 1200-1300 K and 4.0 dex for the T component. This is dif¬ 
fering temperatures and gravities, which agrees with the results 
of lFahertv et al.l(l2014l) . As previously discussed, given the simi¬ 
lar luminosities of the two pair members and assuming coevality, 
the TeffS of the two objects cannot differ by more than ~60 K (at 
the 2 cr level). In the following, we intend to investigate whether 
one single temperature and gravity can provide reasonable fits to 
the observed data. 

The observed differences between the X-shooter and com¬ 
puted BT-Settl spectra show the lack of our understanding of 
the physics of atmospheres of L and T dwarfs in general. To 
improve the fits with a special attention to the simultaneous fit 
of the sodium and potassium subordinate and resonance lines, 
we computed synthetic spectra of L and T dwarfs following 
the procedure by iPavlenko et al.l (12007al) . as summarised be¬ 
low. We took advantage of the versatility of the WITA code 
to modify the relative abundances and number densities of 
species shaping the spectral energy distributions of L and T 
dwarfs. We computed chemical equilibrium for a mix of approx¬ 
imately 100 molecul ar spec ies and adopted the conventional ap¬ 
proach of iPavlenko] (12013l) to treat the depletion of atomic al¬ 
kalis into molecular species. Then, we computed partial pres- 
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sures of some molecules and atoms which exceeded the pres¬ 
sures of the gas-dust phase transitions and decreased their gas 
phase abundances to the corresponding equilibrium values (see 
iPavlenkoii 19981) . We computed the profiles of the Nal and Kl 
resonance d oublets in the framewor k of a quasi-static approach 
described in Pavlenko et al.l ( 2007bl ) with an upgraded approach 
from lBurrows et ah ( 20031). We computed the CrH and FeH line 
lists fr om the work of Burrows et al.l (l2002bl) and iDulick et al.l 
(120031). respectively, while t he water vapour line list was taken 
from lTennysonet aD d2007h . We considered the spectroscopic 
data for atomic absorption of the Vienna Atomic Line Database 
(VALD; iKupka et al.l 11999iP l We plot the overall spectral en¬ 
ergy distribution produced by WITA after modifying the num¬ 
ber densities of Nal, Kl, CrH, FeH, and water vapour for L and 
T component in the bottom plots in Fig. [9] assuming an effec¬ 
tive temperature of HOOK and a gravity of log(g) of 5.0-5.5. 
In that way we use rather phenomenological approach to fit our 
theoretical spectra to the observed spectral energy distributions. 
The detailed analysis of these abundances deserve a more dedi¬ 
cated study which is beyond the scope of this paper. Nonetheless, 
it emerges that the spectral properties of L/T transition objects 
are not only determined by effective temperature and gravity but 
also depend on the depletion of potassium and sodium atoms 
into condensates and dust species, as previously suggested by 
iMarlev et al.l (i2010h . 


7. Summary 

We presented high-quality optical and near-infrared of the 
Luhman 16AB obtained with VLT/X-shooter. We classify each 
component of this nearby brown dwarf binary as L6-L7.5 and 
T0+1, respectively. We focussed our analysis primarily on the 
strength of alkali lines, which can be summarised as follows: 

• We measured similar pEWs for the lithium absorption line 
in each component, placing both objects in the substellar 
regime. The pEW measurement of the L component is con¬ 
sistent with the values of other late-L field dwarfs. Lithium is 
detected for the first time in a T dwarf (see also lFahertv et alJ 
l20l4h with an pEW of 8.4+1.5A 

• We set an upper limit of 30 mA on the lithium subordinate 
line at 812.64 nm. 

• The pEWs of the rubidium lines are similar for the L and T 
components of the binary whereas the cesium lines appear 
stronger in the T dwarf. 

• We detected the sodium subordinate lines in the red optical 
spectra of the L and T components. The L-type member ex¬ 
hibits stronger sodium lines than the T component. 

• The pEWs of lithium and cesium decrease for spectral types 
later than ~L6 and T0-T5, respectively, likely correspond¬ 
ing to the condensation of these species into molecules in 
the atmospheres of cool brown dwarfs. Our comparison with 
other L and T dwarfs shows that the condensation of cesium 
occurs at temperatures about 300 K cooler than the one of 
lithium, as predicted by models. 

We derived effective temperatures of 1305^K and 1320 + K 
and luminosities of -4.68 + 0.08 dex and -4.66 ±0.13 dex for 
the L and T component of the binary, respectively, which are 
in agreement with the literature. Moreover, we measured helio¬ 
centric radial velocities of 25.1 and 19.3 kms~' for each com¬ 
ponent, yielding Galactic space velocities which do not indicate 

5 http ://vald .astro, uni vie. ac. at/ ~ vald/php/vald. php 


membership to any of the known, nearby, young moving group. 
Finally, we fit synthetic spectra to the observed spectra of each 
component of the binary by combining the BT-Settl models and 
the WITA chemical abundance code. We reproduced satisfacto¬ 
rily the spectral energy distributions of the L and T dwarfs with 
a single set of temperature and gravity, but different relative at¬ 
mospheric abundances of some elements like K, Na, Cr, which 
may be related to varying levels of dust condensation. 
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Fig. 9. Top: Comparison of the X-shooter optical and near-infrared spectra (black line) of the L (left) and T (right) components of 
Luhman 16AB and the BT-Settl model atmosphere with 1400/5.0 (red line) and normalised around 800 nm to match the potassium 
wings. Bottom: Same X-shooter spectra (black line) compared to synthetic spectra generated combining the WITA chemical code 
with the BT-Settl model atmosphere for a temperature of 1400 K, a gravity of log(g) of 5.0-5.5, and different relative abundances of 
Na, K, and other refractory elements (smoothed for visualisation purposes). 
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